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STRESS ANALYSIS OF LAMINATED GRAPHITE/EPOXY COMPOSITE PLATE USING 



This work presents a stress analysis of Graphite/Epoxy laminated composite plate. In the present work the stress 
behaviour of laminated composite plates under Tensile loading using a four-node element with six degrees of freedom at 
each node and translations in the x and y directions. The static stress analysis includes the all type of stress behaviour in 
diagrammatic form and results are plotted for investigation. In the present study the modelling is done in ANSYS 14.5 .In 
this study investigations were carried on square plates starting with Three layers of top location of 0° angle ply laminated 
composite plates at clamped boundary condition and 0° angle ply oriented with 15,30, 45,60,and 75 degree then analyze 
the stresses for optimizing the structure. 
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1. INTRODUCTION 

Composite materials are extensively used in aerospace, automobile, nuclear, marine, biomedical and civil 
Engineering. As composite materials having high strength to weight ratio and high stiffness to weight ratio. The material 
has superior fatigue characteristics and also ability to change fibre orientations to meet design requirements. Now a day's 
laminated composite materials are the primary need of high rise buildings. In practical the use of laminated composite as 
slabs, beams, panels, deck etc. The matrix, compared with fibre direction, limits the strength of laminated composites. 
However composite structures subjected to low- velocity impacts or the drop of minor objects, such as tools during 
assembly or maintenance operation, exhibit a brittle behaviour and can sustain significant damage. 

These impacts are particularly significant in design of laminated plies. The structural components like beam made 
of composite materials are being increasingly used in engineering applications. Because of their complex behaviour in the 
analysis of such structures some technical aspects must be taken into consideration. Finite element method is versatile and 
efficient tool for the analysis of complex structural behaviour of the laminated composite structures. The analysis of 
vibration and dynamics, buckling and post buckling, failure and damage analysis based on the various laminated plate 
theories is mainly carried out using ANSYS software. 

2. LITERATURE REVIEW 

Wimmer and Pettermann [1] numerically simulated separation in composite plates. In his study, he used 
combination of critical and failure loadings. Yang et al. [2] defined exact numerical 
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method for mechanical behaviour of composite plates with elastic joints. They assumed isotropic joints and also considered 
friction in contact surface between joints and plates. In their analysis, plates were assumed symmetric. L. Jun, H.Hongxing 
and S. Rongying [3], introduced a dynamic finite element method for free vibration analysis of generally laminated 
composite beams on the basis of first-order shear deformation theory. The influences of Poisson effect, couplings among 
extensional, bending and tensional deformations, shear deformation and rotary inertia are incorporated in the formulation. 
The dynamic stiffness matrix is formulated based on the exact solutions of the differential equations of motion governing 
the free vibration of generally laminated composite beam. Hyer and Liu [4], investigated effect of pin elasticity, clearance 
and friction on radial and tangential stress distribution around hole in mechanical joints of orthotropic plates by numerical 
method. Their model was 2D and constituted of two symmetric layers to investigate aforementioned parameters. Yegao et 
al [5] presented a general formulation for free and transient vibration analysis of composite laminated beams with arbitrary 
lay- ups and any boundary conditions. A modified variational principle combined with a multi-segment partitioning 
technique is employed to derive the formulation based on a general higher order shear deformation theory. The material 
coupling for bending- stretching, bending-twist, and stretching twist as well as the poison's effect are taken into account. 
K. Mosalam, L Glascoe, J Bernier [6] stated as different varieties of bricks such as hollow bricks, fly ash bricks, clay 
bricks, wire cut bricks, etc. and varied mortar mixes containing fly ash, lime, or only cement and sand. Other factors 
affecting the properties of masonry are anisotropy of the brick units, mortar joint thickness, mortar strength and 
workmanship. Therefore, the masonry properties are not similar and vary with the properties of its constituent materials 
like brick unit and mortar. To understand the behaviour of masonry, the material properties of brick units and mortar has to 
be determined. Shahid Nazir, Manicka Dhanasekar, [7] has modelled the masonry for thin layered mortar joints. 
The reviews helped in carryout the present work. 

3. METHODOLOGY 

3.1 Finite Element Method (FEM) 

The physical structure that was used in this work is a fibre reinforced composite plate; the length (a) of 1.0 metre 
width (b) of 0.5metre and thickness (h) of the plate is 0.015m. A number of analyses are performed in this design study, 
using a finite element model of the plate. The model was developed using linear layered structural shell elements in 
ANSYS 14.5. The global x-coordinate is taken along the length of the plate; the global y-coordinate is taken along the 
width of the plate while the global z-direction is taken out the plate surface. There are 40 elements in the axial direction 
and 40 along the width one. The boundary conditions, all sides are constrained in all directions. The pressure applied on 
the plate is 1000 N. In this study, three ply [0°/30 0 /-45°] symmetric laminated composite plate is considered in the analysis. 
The plate is analyzed for deflections and stresses under a fully constrained boundary condition when the plate is subjected 
to tensile loading along the X and Y - directions for various ply locations. 

3.2 Geometry of the Shell Element 

In ANSYS software, there are many element types available to model layered composite materials. In our FE 
analysis, the linear layered structural shell element is used. It is designed to model thin to moderately thick plate. An 
accurate representation of irregular domains (i.e. domains with curved boundaries) can be accomplished by the use of 
refined meshes and/or irregularly shaped elements. For example, a non-rectangular region cannot be represented using only 
rectangular elements; however, it can be represented by triangular and quadrilateral elements. Since, it is easy to derive the 
interpolation functions for a rectangular element and it is much easier to evaluate the integrals over rectangular geometries 



Impact Factor (JCC): 3.6234 



NAAS Rating: 2.02 



Stress Analysis of Laminated Graphite/Epoxy Composite Plate Using Fem 



11 



than over irregular geometries, it is practical to use quadrilateral elements with straight or curved side assuming a means to 
generate interpolation functions and evaluate their integrals over the quadrilateral elements. The linear layered structural 
shell element is shown in Figure 3. 1 Nodes are represented by I, J, K, L, M, N, O, and P. 




Figure 1: Geometry of 8-Node Element with Six Degrees of Freedom 



In the present context of table (3.1&3.2), modelling the stresses in a Symmetric Cross-Ply Laminar Composite in 
tension load, in ANSYS Mechanical APDL. This composite consists of three layers. The model will use two dimensional 
layered shell elements. By comparing the results with the analytical solution based on the Generalized Hooke's Law, the 
module will emphasize techniques on modelling orthotropic materials and layered materials. 

3.3 Problem Description 

An orthotropic plate with three numbers of layers is subject to transverse loading condition for clamped boundary 
condition and the best ply location oriented with different angles has been considered for the present study, and the results 
were given in diagrammatic form. 



Table 3.1: Geometric Properties of Orthotropic Material Plates 



Dimensions 


No. 
of Layers 


Stacking 
Sequence 


Fiber Orientation 


Length=1.0m 

Width=1.0m 
Thickness=0.015m 


3 


0°/30 0 /-45° 


Symmetric ply 



Table 3.2: Material Properties of Graphite/Epoxy Composite Material 



Elastic Constants 


Values 


Ei 


181GPa 


E 2 


10.3 GPa 


E 3 


10.3 GPa 


v 12 =v 13 


0.28 


v 23 


0.01593 


Gi2=Gi3 


7.17 GPa 


G 2 3 


7.17 GPa 
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RESULTS & DISCUSSIONS 

The problem was anlyzed using FEM based software ANSYS. Figures detail the different stresses of laminated 
composite plate. The stress in X direction, stress in Y direction, Shear stress in XY direction for 0° layer located at bottom, 
middle and top location of the lamina are shown in Figure 4.1 to 4.9. The orientation of layers are changed and change 
values are plotted in Figure 4. 10 to 4.24. The optimized ply location and ply orientation are identified with the results. 



(a) When Layer is Located at Bottom (0 ) 



(b.) When Layer is Located at Middle(0 u ) 




Figure 4.1: Stress in X Direction 



Figure 4.4: Stress in X Direction 




Figure 4.2: Stress in Y Direction 



Figure 4.5: Stress in Y Direction 




Figure 4.3: Shear Stress in XY Direction 



Figure 4.6: Shear Stress in XY Direction 



(d) When 0° Layer is Oriented with 15° 



(c) When Layer is Located at Top (0 U ) 




Figure 4.7: Stress in X Direction 



Figure 4.10: Stress in X Direction 
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Figure 4.8: Stress in Y Direction 




Figure 4.11: Stress in Y Direction 



Figure 4.9: Shear Stress in XY Direction 



Figure 4.12: Shear Stress in XY Direction 



(e) When 0° Layer is Oriented with 30° 



(f) When 0 U Layer is Oriented with 45° 





Figure 4.13: Stress in X Direction 



Figure 4.16: Stress in X Direction 





Figure 4.14: Stress in Y Direction 



Figure 4.17: Stress in Y Direction 




Figure 4.15: Shear Stress in XY Direction 



Figure 4.18: Shear Stress in XY Direction 
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(g) When 0° Layer is Oriented with 60° (h)When 0° Layer is orienEted with 75 




Figure 4.19: Stress in X Direction Figure 4.22: Stress in X Direction 




Figure 4.20: Stress in Y Direction Figure 4.23: Stress in Y Direction 




Figure 4.21: Shear Stress in XY Direction Figure 4.24: Shear Stress in XY Direction 



Graph 4.1 Stresses Vs Ply Locations 



STRESSES Vs PLY LOCATIONS 




Stress in Y direction 
Shear Stress in XY direction 



Top Middle Bottom -^-Stress in X direction 

PLY LOCATION 



Graph 4.2 Stresses Vs Ply Orientations 
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STRESSES Vs PLY ORIENTATION 




15 30 45 60 75 

PLY ORIENTATION (Degrees) 



Shear Stress in XY Direction 
•Stress in Y Direction 
•Stress in X Direction 



Graph 4.3 Strains Vs Ply Locations 



Table 4.1: Tabulation of Results: Regarding the Location of 0° Ply 



Design Parameter 


When 0° Layer is 
Located at Bottom 


When 0° Layer is 
Located at Middle 


When 0° Layer is 
Located at Top 


Stress in X Direction(Pa) 


5.5767xl0 4 


lO.lOlxlO 4 


3.3513xl0 4 


Stress in Y Direction(Pa) 


4.5312xl0 4 


5.0510xl0 4 


6.1875xl0 4 


Shear Stress in XY 
Direction (Pa) 


-1.2800xl0 4 


-9.067X10 4 


-2.7504xl0 4 



Different ply location and respective stresses are detailed in Table 4.1, The bottom ply location is the optimized 
location for stresses in Y direction and for the shear stress. 

4.2 Optimization of Orientation Ply Orientation 

The 0° location at bottom ply location is oriented in 15 0 ,30 0 ,45 0 ,60°,75 0 and the values of stresses are developed 
to estimate optimised ply orientation in tensile loading of a laminated composite plate and the results are tabulated in 
Table 4.2. 

Table 4.2 Tabulation of Results Regarding Optimization of Ply Orientation 



Design Parameter 


When 0° 
Layer is 
Oriented 
with 15° 


When 0° Layer 
is Oriented 
with 30° 


When 0° Layer 
is Oriented with 
45° 


When 0° 
Layer is 
Oriented 
with 60° 


When 0° Layer is 
Oriented with 

75° 


Stress in X 
Direction (Pa) 


1.7823xl0 4 


6.9241 xlO 3 


7.6931xl0 3 


1.8365xl0 4 


2.9591xl0 4 


Stress in Y 
Direction (Pa) 


5.8004xl0 4 


4.8094xl0 4 


2.9193xl0 4 


7.537xl0 3 


-5.230xl0 3 


Shear Stress in XY 
Direction(Pa) 


-3.0518xl0 4 


-4.0815xl0 4 


-5.1469xl0 4 


5.3214xl0 4 


-4.4121xl0 4 



5. CONCLUSIONS 

In this paper, Finite element analyses have been carried out for various ply locations (0°, 30°,- 45°) and 
orientations. The study the effect of stresses of laminated composite plates subjected to tensile load. From the results it is 
concluded that The stress in Y direction and shear stresses are minimum if ply location is at bottom. Also concluded for 
load application in X direction 30° is optimum orientation of the composite Lamina. The stresses Y direction is minimum 
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at 45°and for shear stress is 15° for the same loading conditions indicates the load carrying property changes with the ply 
orientation of the composite. 

5.1 Future Scope 

The work may be extended by changing the aspect ratio, length to width ratio (a/b), width to thickness ratio(b/h) 

also.. 
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